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Q\ ' Abstract 

■ Extensive statistical-equilibrium calculations on neutral sulfur and zinc were car- 

O ! ried out, in order to investigate how the non-LTE effect plays a role in the determina- 

^ \ tion of S and Zn abundances in F, G, and K stars. Having checked on the spectra of 

representative F-type stars (Polaris, Procyon, and a Per) and the Sun that our non- 
LTE corrections yield a reasonable consistency between the abundances derived from 
different lines, we tried an extensive non-LTE reanalysis of published equivalent-width 
data of S I and Zn I lines for metal-poor halo/disk stars. According to our calculations, 
S I 9212/9228/9237 lines suffer significant negative non-LTE corrections amounting to 
< 0.2-0.3 dex, while LTE is practically valid for S I 8683/8694 lines. Embarrassingly, 
as far as the very metal-poor regime is concerned, a marked discordance is observed 
between the [S/Fe] values from these two abundance indicators, in the sense that 
the former attains a nearly flat plateau (or even a slight downward bending) while 
the latter shows an ever-increasing trend with a further lowering of metallicity. The 
reason for this discrepancy is yet to be clarified. Regarding Zn, we almost confirmed 
the characteristic tendencies of [Zn/Fe] reported from recent LTE studies (i.e., an 
evident /slight increase of [Zn/Fe] with a decrease of [Fe/H] for very metal-poor/disk 
stars), since the non-LTE corrections for the Zn I 4722/4810 and 6362 lines (tending 
to be positive and gradually increasing towards lower [Fe/H]) are quantitatively of 
less significance (^$0.1 dex). 

Key words: Galaxy: evolution — line: formation — stars: abundances — stars: 
atmospheres — stars: late-type 
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1. Introduction 

The subject of this paper is to investigate how the non-LTE effect influences the spec- 
troscopic determination of sulfur and zinc abundances in F-, G-, and K-type stars (especially 
metal-poor ones) used for studying the Galactic chemical evolution history. 

1.1. Astrophysical Significance of S and Zn 

Sulfur belongs to the group of "a-capture elements" (along with O, Mg, Si, Ca, and 
Ti), a large fraction of which are considered to have been synthesized in short-lived massive 
stars (and thrown out by type II supernovae) at the early-time of the Galaxy. On the other 
hand, while the production of zinc is still controversial and several possibilities are discussed 
(see, e.g., a summary presented by Chen et al. 2004), it is often regarded as being associated 
with "Fe-group elements" which were significantly produced in the later stage of the Galactic 
history by type la supernovae of longer-lived intermediate-mass stars. 

The particularly important characteristic of these two elements is that they are chem- 
ically "volatile" such as the cases of C, N, and O. That is, owing to their low condensation 
temperature (T c ~ 650 K), they are difficult to condense into solids, unlike other "refractory" 
elements (Mg, Si, Fe-group elements, etc.) with a high T c of ~ 1300-1500 K. Namely, S and Zn 
are considered to be hardly affected by depletion due to dust formation. This fact is especially 
significant in the chemical composition of intergalactic gas, where these two volatile elements 
are likely to retain their original composition even for such a condition, while other refractory 
species (such as Mg or Fe) may have been significantly fractionated onto dust and depleted. 
For this reason, in the analysis of damped Lyman a (DLA) system of QSO absorption lines, 
S and Zn are generally regarded as being (depletion-independent) important tracers of the a 
group and the Fe group, respectively, which provides us a possibility to use [S/Zn] and [Zn/H] 
determined from DLA as a "chemical clock" of high- z universe (see, e.g., the summary of Nissen 
et al. 2004 and the references therein). Anyway, as a first step toward such an advanced appli- 
cation, the nucleosynthesis history of these elements in our Galaxy has to be well understood 
by observationally establishing the behavior of [S/Fe] and [Zn/Fe] with a change of [Fe/H] in 
metal-poor stars. Yet, this problem has not necessarily been straightforward. Especially, we 
are still in a confusing situation concerning the case of S. 

1.2. Controversy over the Behavior of [S/Fe] 

Since sufficiently strong sulfur lines (measurable even in late-type metal-poor stars) are 
located in the near-IR region, earnest investigations on the [S/Fe] vs. [Fe/H] relation began 
in the 1980's when efficient solid-state detectors had become widely used. Further, these early 

* The electronic tables (El, E2, E3, and E4) will be made available at the E-PASJ 
web site upon publication, while they are provisionally accessible at the WWW site of 
(http://optik2.mtk.nao.ac.jp/~takeda/sznabund/). 
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studies based on S I 8693/8694 lines sufficed to reveal the gradual increase of [S/Fe] from ~ 

(at [Fe/H] ~ 0) to ~ +0.5 (at [Fe/H] 1.5) with a decrease of metallicity (Clegg et al. 1981; 

Frangois 1987, 1988), which is typical for the a group. 

The more important issue is, however, how it behaves itself in the very metal-poor regime 
of —3 < [Fe/H] < —2. The first two trials of investigating [S/Fe] in this region based on the S I 
8683/8684 lines (being very weak) resulted in a conclusion that [S/Fe] continues to ever increase 
with a further lowering of [Fe/H] (up to [S/Fe] ~ +1 at [Fe/H] ~ —3), while retaining its slope 
almost unchanged (Israelian, Rebolo 2001; Takada-Hidai et al. 2002). This makes a marked 
contrast with other a-group elements (e.g., Mg, Si, Ca, Ti) showing a "flat" tendency (i.e., 
[a/Fe] attains an almost constant value below [Fe/H] < — 1). Thus, these authors suggested 
that [S/Fe] (along with [O/Fe] 1 ) does not conform to the standard a behavior, and that an 
adequate chemical evolution model (e.g., hypernovae or time-delayed deposition) reasonably 
reproducing such a tendency may be required. 

Then, the situation turned into a new aspect, when the S I lines of the 9212/9228/9237 
triplet (much stronger than S I 8693/8694 lines) began to be used for exploring the sulfur 
abundances of very metal- deficient stars based on the two recent studies of Ryde and Lambert 
(2004) and Nissen et al. (2004). Interestingly, they both arrived at the same conclusion that 
[S/Fe] attains an approximately constant value of ~ +0.3 at [Fe/H] < — 1, which is nothing 
but a typical a-group trend explainable with a standard chemical evolution model. Namely, 
their results derived from the 9212/9228/9237 lines are evidently opposed to what Israelian 
and Rebolo (2001) and Takada-Hidai et al. (2002) reported based on the 8693/8694 lines. The 
reason for this discrepancy should be clarified by all means. 

Here, one point should be mentioned. Their analyses are based on the belief that reliable 
abundances can be obtained from the 9212/9228/9237 lines on the assumption of LTE. It 
appears, however, that their arguments are a kind of speculation, lacking any firm background; 
namely, they seem to overinterpret the consequence of Takada-Hidai et al. (2002) (i.e., non- 
LTE corrections are negligible for the 8693/8694 lines of multiplet 6) and extend it to the 
9212/9228/9237 lines of multiplet 1, without paying attention to the large difference in the 
line-formation mechanism between these two multiplets. Therefore, in order to examine the 
validity of their treatment, an extensive study has to be carried out on how and to what extent 
the non-LTE effect is important for these multiplet 1 sulfur lines over a wide range of stellar 
parameters. 

In addition, a comprehensive comparison study for the lines of both multiplets (1 and 6) 
on a data sample comprising as many stars as possible is necessary. Takada-Hidai et al. (2005) 
very recently published their new observations of both 8693/8694 and 9212/9228/9237 lines 
for some ~ 20 stars carried out at Okayama Astrophysical Observatory, and made an effort to 

1 Oxygen is also suspected to show such an ever-increasing trend though it still remains as a controversial 
matter (i.e., different behaviors are suggested from different lines; cf. Takeda 2003 and the references therein). 
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investigate how the abundances derived from these two different multiplets compare with each 
other, while invoking the assumption of LTE as was done by Ryde and Lambert (2004) and 
Nissen et al. (2004). Unfortunately, however, the insufficient number of data at [Fe/H] < —2 
in their sample, coupled with a large scatter, did not allow any convincing argument on the 
[S/Fe] vs. [Fe/H] relationship of metal-poor stars or on the validity of LTE (cf. their figure 3). 
This lesson actually taught us the necessity of carrying out an extensive non-LTE reanalysis on 
almost all available observational data published so far combined together, which is thus the 
first purpose of this investigation. 

1.3. Does Zn Scale with Fe? 

The run of [Zn/Fe] with [Fe/H] in Galactic metal-poor stars is not as controversial as 
the case of sulfur. Nevertheless, some problematic points have emerged from recent studies. 

Since the pioneering work of Sneden and Crocker (1988) and Sneden, Gratton, and 
Crocker (1991), it has long been believed that Zn practically scales with Fe (i.e., [Zn/Fe] ~ 0) 
over a wide range of metallicities (—3 ^ [Fe/H] < 0). As a matter of fact, this is the reason 
why Zn is regarded as a substitute for Fe in the condition where the refractory Fe may be 
appreciably depleted (e.g., analysis of DLA line system mentioned in subsection 1.1). 

However, recent studies extending to an ultra-metal-poor region, or those with higher 
accuracy than before, have gradually cast doubt on this simple picture. For example, at metal- 
licities lower than [Fe/H] < —3, it has been reported that [Zn/Fe] progressively increases with a 
decrease in [Fe/H] (cf. Primas et al. 2000; Cayrel et al. 2004). Also, according to Chen, Nissen, 
and Zhao (2004) investigation on the Zn abundances of Galactic disk stars based on the weaker 
Zn I 6362 line, which may be more reliable than the usually used 4722/4810 lines, [Zn/Fe] ap- 
pears to show a slightly increasing tendency as [Fe/H] decreases from ~0 to ~ —0.6, consistently 
with the findings of Reddy et al. (2003) and Bensby, Feltzing, and Lundstrom (2003). Bensby, 
Feltzing, and Lundstrom (2003) also suggested a tendency of "uprising" [Zn/Fe] at [Fe/H] > 
as well as a separation of [Zn/Fe] between thick-disk and thin-disk stars. Furthermore, the 
[Zn/Fe] vs. [Fe/H] relation constructed from published results (combined) appear to suggest 
a shallow/marginal dip at [Fe/H] ~ —1.2 (cf. figure 8 of Mishenina e al. 2002 or figure 7 of 
Nissen et al. 2004; though the number of stars is rather insufficient around this metallicity.) In 
any case, there are now several good reasons to suspect that the real situation is not so simple 
(i.e., scaling with Fe) as has ever been considered. Hence, it is needed to establish the trend of 
[Zn/Fe] as a function of metallicity as precisely as possible. 

Here, it is one of our concerns that all of these studies were done under the assumption 
of LTE. As a matter of fact, no investigation has so far been carried out on the non-LTE effect 
of zinc, so far as we know. Then, do some of the departures from [Zn/Fe] ~ described above 
have something to do with neglecting the non-LTE effect? Or alternatively, is the non-LTE 
correction so negligible as to guarantee the reliability of such trends derived from LTE? In search 
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of the answers to these questions, we decided to carry out statistical equilibrium calculations on 
neutral zinc, in order to estimate the non-LTE abundance corrections for Zn I 4722/4810/6362 
lines for a wide range of stellar parameters, and to perform an extensive non-LTE reanalysis of 
published equivalent-width data of these Zn I lines toward establishing the [Zn/Fe] vs. [Fe/H] 
relation in our manner. This is the second aim of this study. 

2. Non-LTE Line Formation and Abundance Corrections 

2.1. Sulfur 

2.1.1. Non-LTE calculations on S I atom 

The procedures of our statistical-equilibrium calculations for neutral sulfur are almost 
the same as those described in Takada-Hidai and Takeda (1996) and Takada-Hidai et al. (2002). 
For the present study, however, we reconstructed a new S I model atom comprising 57 terms 
(up to 3s 2 3p 3 8f 3 F at 81837 cm -1 ) and 191 transitions, while using Kurucz and Bell's (1995) 
compilation of atomic data (gf values, levels, etc.), which we believe to be more realistic than 
the previous one (56 terms and 173 transitions) based basically on Kurucz and Peytremann's 
(1975) data. 

The treatment of the photoionization cross sections is the same as described in subsection 
3.3 of Takada-Hidai and Takeda (1996); namely, the available cross-section values compiled by 
Mathisen (1984) were adopted for the lowest three terms (original sources: Tondello 1972 for 
3p 4 3 P; Chapman, Henry 1971 for 3p 4 X D; McGuire 1979 for 3p 4 1 S), while the hydrogenic 
approximation was assumed for the remaining terms. 

Regarding the collision cross section, we followed the recipe adopted in subsubsection 
3.1.3 of Takeda (1991). It should also be mentioned that no correction [i.e., a correction factor 
of k = 1 or logarithmic correction of h(= log A;) = 0] was applied to the H I collision rates 
computed with the classical approximate formula (Steenbock, Holweger 1984; Takeda 1991), 
though test calculations with varying k from 10 to 10~ 3 were also performed (see below). The 
validity of this choice is discussed in connection with the analysis of bright F stars presented 
in section 3. 

Since we planned to make our calculations applicable to stars from metal-rich (population 
I) down to very low metallicity (extreme population II) at early-F through early-K spectral types 
in various evolutionary stages (i.e., dwarfs, subgiants, giants, and supergiants), we carried out 
non-LTE calculations on an extensive grid of 210 (6x5x7) model atmospheres, resulting from 
combinations of six T efr values (4500, 5000, 5500, 6000, 6500, 7000 K), five logg values (1.0, 
2.0, 3.0, 4.0, 5.0), and seven metallicities (represented by [Fe/H]) (+0.5, 0.0, -0.5, -1.0, -2.0, 
—3.0, —4.0). As for the stellar model atmospheres, we adopted Kurucz's (1993) ATLAS9 models 
corresponding to a microturbulent velocity (£) of 2 km s -1 . Regarding the sulfur abundance 
used as an input value in non-LTE calculations, we assumed Ag nput = 7.21 + [Fe/H] + [S/Fe], 
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where the values of 0.0 (for [Fe/H] = +0.5, 0.0, -0.5) and +0.5 (for [Fe/H] = -1.0, -2.0, 
—3.0, —4.0) were assigned to [S/Fe] in order to roughly simulate the behavior of this ratio (cf. 
subsection 4.1). The solar sulfur abundance of 7.21 was adopted from Anders and Grevesse 
(1989) (which is used also in the ATLAS9 models). The microturbulent velocity (appearing in 
the line-opacity calculations along with the abundance) was assumed to be 2 km s _1 , to make 
it consistent with the model atmosphere. 

2.1.2. Non-LTE characteristics of S I line formation 

In figure 1 are shown the S-^{j) / B{t) (the ratio of the line source function to the Planck 
function, and nearly equal to ~ b u /b\, where b\ and b u are the non-LTE departure coefficients for 
the lower and upper levels, respectively) and ^ lte (t)//q TE (t) (the NLTE-to-LTE line-center 
opacity ratio, and nearly equal to ~ b\) for each of the multiplet 1 (9212/9228/9237) and mul- 
tiplet 6 (8693/8694) transitions for a representative set of model atmospheres. Two especially 
important characteristics can be read from this figure: 

- Generally, the inequality relations of S^/B < 1 (dilution of line source function) and 
^nlte^lte > ^ (enhanced line-opacity) hold in the important line-forming region for both cases 
of multiplets 1 and 6, which means that the non-LTE effect almost always acts in the direction 
of strengthening the 9212/9228/9237 and 8693/8694 lines; i.e., the non-LTE correction is gen- 
erally negative. 

- The departure from LTE in the line opacity (the enhancement of ^ lte /Zq TE over 1) becomes 
prominent for higher T c g and/or very low-metallicity case, which indicates that the non-LTE 
effect may become significant in very low metallicity regime of [Fe/H] ~ —3 down to —4 even 
if the line-strength is weak. 

2.1.3. Grid of non-LTE corrections for S I lines 

Based on the results of these calculations, we computed extensive grids of theoretical 
equivalent-widths and the corresponding non-LTE corrections for the eight selected important 
lines (S I 8693 and 8694 lines of multiplet 6; 9212, 9228, and 9237 lines of multiplet 1; 10455, 
10456, and 10459 lines of multiplet 3) for each of the model atmospheres as follows. 

For an assigned sulfur abundance (A a ) and microturbulence (£ a ), we first calculated 
the non- LTE equivalent width (W NLTE ) of the line by using the computed non-LTE departure 
coefficients (b) for each model atmosphere. Next, the LTE (A L ) and NLTE (A N ) abundances 
were computed from this M/ NLTE while regarding it as if being a given observed equivalent 
width. We could then obtain the non-LTE abundance correction, A, which is defined in terms 
of these two abundances as A = A N — A h . 

Strictly speaking, the departure coefficients [b(r)} for a model atmosphere correspond to 
the sulfur abundance and the microturbulence of A™ 1511 * and 2 km s _1 adopted in the non-LTE 
calculations (cf. subsubsection 2.1.1). Nevertheless, considering the fact that the departure 
coefficients (i.e., ratios of NLTE to LTE number populations) are (unlike the population itself) 
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not very sensitive to small changes in atmospheric parameters, we also applied such computed 
b values to evaluating A for slightly different A a and £ a from those fiducial values assumed in 
the statistical equilibrium calculations. Hence, we evaluated A for three A a values (/lg nput and 
±0.3 dex perturbations) as well as three £ values (2 km s _1 and ±1 km s _1 perturbations) for 
a model atmosphere using the same departure coefficients. 

We used the WIDTH9 program (Kurucz 1993) for calculating the equivalent width for 
a given abundance, or inversely evaluating the abundance for an assigned equivalent width. 
Actually, this program was considerably modified in many respects: e.g., the treatment of a 
blended feature due to multiplet components, the incorporation of the non-LTE departure in 
the line source function as well as in the line opacity, etc. The adopted line data (gf values, 
damping constants, etc.) are given in table 1. Since the S I 9228 lines are located on the wing 
of Paschen line (H I P 9 9229.0), we replaced the hydrogen-line opacities of the original Kurucz's 
code [based on classical Griem's (1960, 1967) approximation] by a more updated one based on 
the extended VCS theory calculated by Lemke (1997). 

As a demonstrative example of non-LTE corrections, we give the £ = 2 km s _1 results 
for the S I 8694 and 9212 lines computed for representative parameters in table 2, where we 
also present the cases of h = +1,-1,-2,-3 in addition to the fiducial h = [h (= log A;) is 
the logarithmic H I collision correction to be applied to the classical formula; cf. subsubsec- 
tion 2.1.1) for comparison. As can be seen from table 2, the non-LTE effect becomes more 
appreciable with a decrease of h, as expected. Also, conspicuously large (negative) non-LTE 
corrections (accompanied by large line-strengths) seen in low-gravity and/or high-T fj stars are 
worth noting. 

Since the S I photoionization cross sections that we adopted may not be sufficiently up 
to date (cf. subsubsection 2.1.1), we investigated how changing the cross-section values by 
factors of 0.1 and 10 would affect the non-LTE corrections; the results are also given in table 
2 ( 5(A)- and 5(A) + ). While the resulting changes are not necessarily straightforward, we 
can see that the extent of the negative corrections tends to be reduced (i.e., less negative) by 
increasing the photoionization, which may be interpreted as that lines tend to be weakened (i.e., 
bringing the non-LTE correction in the positive direction) by the photo-overionization effect. 
From a quantitative view, however, the variations in A are not significant in most cases. Yet, 
the exceptional cases are high-T e ff as well as low- log g stars, for which changes are appreciably 
large to be several-tenth dex (reflecting the importance of photoionization by UV radiation). 

The complete results of the non-LTE corrections (for all combinations of T e g, log g, and 
£ values for each of the 8 S I lines, though only for the case of h = 0) are given in electronic 
table El. 
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2.2. Zinc 



2.2.1. Non-LTE calculations on Zn I atom 

Again invoking Kurucz and Bell's (1995) compilation of atomic data, we constructed a 
Zn I model atom consisting of 44 terms (up to 3d 10 4s 14d 3 D at 75112 cm -1 ) and 87 transi- 
tions. The hydrogenic approximation was assumed for the photoionization rates from all terms. 
Regarding the collisional rates (due to electron and neutral hydrogen), we followed the classical 
formulae described in subsubsection 3.1.3 of Takeda (1991). Although we eventually adopted 
h = (i.e., without applying any correction to the classical value) for the H I collision rates, 
the effect of varying h was also examined (cf. table 3). 

Similarly to the case of S I, we carried out extensive non-LTE calculations on a grid 
of 210 model atmospheres. The input Zn abundance in non-LTE calculations was assumed to 
be A^ ut = 4.60 + [Fe/H] + [Zn/Fe], where we assigned the values of 0.0 (for [Fe/H] = +0.5, 
0.0, -0.5, -1.0, -2.0) and +0.5 (for [Fe/H] = -3.0, -4.0) to [Zn/Fe] while considering the 
recently observed supersolar ratio at extremely low metallicities (cf. subsection 4.2). The solar 
zinc abundance of 4.60 was adopted from Anders and Grevesse (1989). 

2.2.2. Non-LTE characteristics of Zn I line formation 

Similarly to figure 1, we show in figure 2 the behaviors of Si,(t)/B(t) and 
lo LTE (r)/lo TE (r) for two important Zn I transitions of multiplet 2 (4722/4810) and multiplet 6 
(6362) for representative model atmospheres. The noteworthy characteristics recognized from 
this figure are as follows: 

- Unlike the case of S I, S*l tends to be superthermal (Sl > B) at important line-forming 
regions, which generally acts in the direction of line-weakening. 

- The behavior of the NLTE-to-LTE line-opacity ratio differs from case to case; it tends to 
be greater than unity (line-strengthening) at lower T e fj, while it becomes appreciably less than 
unity (line- weakening) at higher T c g (especially for very metal-poor cases). Also, the trend of 
multiplet 2 is significantly different from the case of multiplet 6. 

- Combining these characteristics mentioned above, we can expect that the net non-LTE effect 
is rather complicated, since two mechanisms may occasionally act in the opposite direction and 
compensate with each other. Yet, we may roughly state that the non-LTE effect on the Zn I 
lines is in many cases a slight line- weakening (i.e., positive non-LTE correction), though it may 
sometimes act to strengthen the line case by case. 

2.2.3. Grid of non-LTE corrections for Zn I lines 

As was described in subsubsection 2.1.3 for the case of sulfur, we computed extensive 
grids of the theoretical equivalent-widths and the corresponding non-LTE corrections for the 
three important zinc lines (Zn I 4722 and 4810 lines of multiplet 2; 6362 line of multiplet 6) for 
each of the model atmospheres. The adopted line data are given in table 1. 
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As demonstrative examples of non-LTE corrections, we give the £ = 2 km s _1 results 
for the Zn I 4810 and 6362 lines computed for representative parameters in table 3, where we 
also present the cases of h = +1,-1, —2, —3 in addition to the adopted case of h = 0. Again, it 
is apparent that the non-LTE effect becomes more appreciable with a decrease of h. We see, 
however, the extent of non-LTE correction is generally small and comparatively insignificant; 
also, its sign becomes positive as well as negative case by case. This can be understood from the 
characteristics of non-LTE line-formation described in subsubsection 2.2.2. Yet, appreciably 
large positive non-LTE corrections amounting up to ~0.3 dex seen in low-gravity/high-T c fj/low- 
[Fe/H] cases are worth noting. 

As was done for sulfur, we also investigated how changing the Zn I photoionization cross 
sections affects the non-LTE corrections; the results are given table 3. While we can observe a 
roughly similar tendency to the case of S, the variations are quantitatively insignificant 0.1 
dex). 

The complete results of the non-LTE corrections (for all combinations of T e g, logg, £ 
values for each of the three Zn I lines, though only for the fiducial case of h — 0) are given in 
electronic table E2. 

3. Analysis of S and Zn Lines for Selected F-type Stars and the Sun 

3.1. Importance of Abundance Consistency Check 

According to the calculations described in the previous section, unlike the case of zinc 
where the non-LTE effects on Zn I lines are comparatively minor and insignificant, the extents 
of the non-LTE corrections for S I 9212/9228/9237 lines are appreciably large and of practical 
importance (cf. table 2 in subsubsection 2.1.3 and electronic table El). Thus, it is requisite 
to examine the validity of our calculations (e.g., whether our choice of h = is reasonable 
or not) by empirically checking the abundance consistency between the S I 8693/8694 and 
9212/9228/9237 lines in actual stars. 

Towards this aim, we decided to carry out abundance analyses with respect to these S I 
lines based on our non-LTE calculations, while using the spectra of representative bright F-type 
stars: a Per (F5 lb), Polaris (F8 lb), and Procyon (F5 IV) [along with the solar flux spectra 
of Kurucz et al. (1984)]. The reason why we selected these F stars with T e g > 6000 K is that 
such stars (especially low-gravity supergiants) may serve as the most suitable touchstone for 
this purpose, since the non-LTE effect on these S I lines becomes so large that their sensitivity 
to changing h may be quantitatively appreciable in such a condition of higher T e s and lower 
logg (i.e., the case where lines are strong). 

In addition, we also analyzed the Zn I 4722/4810 and 6362 lines of these F stars and 
the Sun, though their non-LTE corrections are so small that we can not say much about the 
adequacy of our non-LTE calculations. 
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3.2. Observational Data 

The observations of these three target stars (a Per, Polaris, and Procyon) were car- 
ried out by using the new high-dispersion echelle spectrograph GAOES (Gunma Astronomical 
Observatory Echelle Spectrograph), which was recently installed at the Nasmyth Focus of the 
1.5 m reflector of the Gunma Astronomical Observatory and can obtain spectra of high wave- 
length resolution (R ~ 70000 for the standard slit width of 1") along with a wide wavelength 
coverage (~ 1800 A by using the 2Kx4K CCD). [See Hashimoto et al. (2002, 2005) for more 
information.] 

For each star, we obtained spectra at three wavelength regions (region G: 4600-6400 A; 
region R: 5900-7600 A; region I: 7600-9350 A). Most of the data were obtained in the observing 
period of 2004 December 14-17, except for the region G (2004 August 31) and region I (2004 
October 27) spectra of a Per. 

The data reduction (bias subtraction, flat-fielding, aperture-determination, scattered- 
light subtraction, spectrum extraction, wavelength calibration, and continuum-normalization) 
was performed using the "echelle" package of IRAF. 2 

Since the spectrum portion including S I 9212/9228/9237 lines contains numerous H 2 
lines originating from Earth's atmosphere, it was divided by the spectrum of 7 Cas (rapid 
rotator) by using the IRAF task "telluric" to remove these telluric lines, which turned out to 
be reasonably successful in most cases. The final spectra of three stars (along with the solar 
flux spectrum for comparison) at the wavelengths corresponding to the relevant S I and Zn I 
lines are shown in figure 3. 

Based on these spectra, the equivalent widths (EWs) of the lines of interest were mea- 
sured by using the software SPSHOW (in the SPTOOL 3 package developed by Y. Takeda) with 
the Gaussian fitting method or the direct-integration method depending on the cases. We did 
not use the S I 9228 line for the three F stars, because it is blended with the strong Paschen 
line (H I Pg) and less reliable (cf. figure 3). Regarding the equivalent width data of the Sun, 
those of S I 8693/8694 and 9228/9237 lines were taken from Takada-Hidai et al.'s (2005) table 
3 (note that the S I 9212 line could not be measured because of being heavily blended with a 
telluric water vapor line), while the others were newly measured from Kurucz et al.'s (1984) 
solar flux spectrum atlas. The finally resulting EW data used in our analysis are presented in 
table 4. 

2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

3 (http:/ /optik2. mtk.nao.ac.jp/~takcda/sptool/) 



10 



3.3. Abundance Results 

The atmospheric parameters of a Per (T e g = 6250 K, logg = 0.90, [X] = [Fe/H] = 0.0, 
£ = 4.5 km s" 1 ) and Polaris (T cfr = 6000 K, logg = 1.50, [X] = [Fe/H] = 0.0, £ = 5.0 km s" 1 ) 
were taken from Takeda and Takada-Hidai (1994), while those of Procyon (T e g = 6600 K, 
logg = 4.00, [X] = [Fe/H] = 0.0, £ = 2.0 km s -1 ) are the rounded values of the original results 
derived by Takeda et al. (2005). Regarding the Sun, we assumed (T c ff = 5780 K, logg = 4.44, 
[X] = [Fe/H] = 0.0, £ = 1.0 km s" 1 ). 

Again, by using the modified WIDTH9 program as in subsubsection 2.1.3, the abun- 
dances of S and Zn were derived from the EW data and the model atmosphere for each star, 
which was constructed from Kurucz's (1993) ATLAS9 model atmospheres grid by interpolat- 
ing in terms of T c g, logg, and [Fe/H]. The resulting NLTE/LTE abundances and the corre- 
sponding NLTE corrections are given in table 4, where roughly estimated values of the mean 
line-formation depth are also presented. 

We should keep in mind that a too-rigorous quantitative discussion is not very meaning- 
ful, especially for a Per and Polaris, because abundance determinations of supergiants involve 
considerable difficulties (large uncertainties in establishing T eff and logg, depth-dependence of 
£, etc.; cf. Takeda, Takada-Hidai 1994). Yet, we can recognize that the large/evident discrep- 
ancies between the LTE abundances of S I 8693/8694 and 9212/9237 (< 1.1 dex for a Per, < 0.7 
dex for Polaris, and < 0.3 dex for Procyon) have been considerably reduced in the case of the 
NLTE abundances (~ 0.2 dex for a Per, ~ 0.1 dex for Polaris 4 , and ~ 0.1 dex for Procyon). 5 
This is surely an encouraging result, which lends support for the practical reliability of our 
non-LTE calculations (we may state, for example, that h = is not a bad choice for the H I 
collision rates of neutral sulfur atoms). 

In view of the fact that the S I 9212/9228/9237 lines of multiplet 1 suffer appreciably 
large non-LTE corrections, we show in figure 4 how the theoretical non-LTE and LTE profiles for 
the representative S I 9237.54 line (computed with the J 4 NLTE value given in table 4) differ from 
each other, and how they match the observed line profile for each star, in order to demonstrate 
the importance of the non-LTE effect. We can see from this figure that the non-LTE line profiles 
satisfactorily reproduce the observations. 

Regarding Zn, the non-LTE effect on the zinc abundance determinations from Zn I 
lines is insignificant 0.1 dex) even in the low-gravity F supergiants (as expected), and thus 
only from the present results alone we cannot comment much on the validity of our non-LTE 
calculations on zinc, though we can argue that our non-LTE abundances (being practically the 

4 As seen from the resulting NLTE abundances, S abundance of Polaris appears to be subsolar by 0.2-0.3 
dex, though we had better be conservative about this, since it might stem from errors (e.g., too high T e g) 
in atmospheric parameters. 

5 Regarding the Sun, the extents of the non-LTE corrections are so small that they can not be used to judge 
the adequacy of our non-LTE calculations. 
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same as the LTE abundances) do not yield any appreciable inconsistency. 

Some remarks on the damping parameters may be due here. The slight differences 
between the LTE solar S abundances of the present study and those of Takada-Hidai et al. 
(2005) are mostly due to the differences in the adopted damping parameters. Namely, in 
contrast to the present treatment, they used the classical formula for the radiation damping, 
and applied the correction of AlogCg = +0.99 to the Cq value computed from the classical 
Unsold's (1955) formula for the van der Waals effect damping, which is equivalent to multiplying 
r^f^ slcal by a factor of 2.5 (a frequently assumed enhancement factor). As a matter of fact, we 
also investigated how the non-LTE abundances would change by multiplying the van der Waals 
damping width (for which we assumed the classical treatment, essentially equivalent to Unsold's 
formula, for all lines; cf. table 1) by a factor of 2.5, as also shown in table 4. As expected, 
appreciable negative variations (amounting to 0.1-0.2 dex) are seen for the high-gravity stars of 
the Sun and Procyon. However, since this increase in T v d w deteriorates the consistency between 
the A NLTE values derived from S I 8693/8694 and 9212/9228/9237, we are reluctant to apply 
such a correction, considering that the classical treatment (as we adopted) is still preferable at 
least for the lines in question. Similarly, almost the same argument holds for the Zn I lines, as 
implied from the consistency between y4 NLTE (4722/4810) and A NL te(6362). 

4. Sulfur Abundances in Metal-Poor Stars 

4-1. Reanaly sis of Literature Data 

We are now ready to study the [S/Fe] vs. [Fe/H] relation of metal-poor stars by re- 
analyzing the published equivalent-width data of the S I 9212/9228/9237 and 8693/8694 lines 
while applying non-LTE corrections based on our calculations. For this purpose, we invoked 
the following papers published so far: Clegg, Lambert, and Tomkin (1981), Frangois (1987, 
1988), Israelian and Rebolo (2001), Takada-Hidai et al. (2002, 2005), Chen et al. (2002), Ryde 
and Lambert (2004), and Nissen et al. (2004). Although our literature survey is not complete, 
we consider that we have picked up most of the important studies, in which the observational 
data are explicitly presented. 

We adopted the same T eff , \ogg, [Fe/H], and £ values as those used in the literature, 
from which the data of the equivalent widths were taken. The Kurucz's (1993) grid of ATLAS9 
model atmospheres and depth-dependent non-LTE departure coefficients were interpolated with 
respect to T eff , logg, and [Fe/H] of each star. Then, as in subsection 3.3, the modified WIDTH9 
program was invoked for determining the non-LTE abundance (A$ LTE ) while using the line data 
given in table 1. Finally, the [S/Fe] ratio was obtained as 

[S/Fe] = (A^ LTE - 7.21) -[Fe/H], (1) 

where 7.21 is the solar sulfur abundance (in the usual scale of log en = 12) taken from Anders 
and Grevesse (1989). 
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In deriving the final [S/Fe] values to be examined, we treated each of the [S/Fe] values 
derived from the S I 9212/9228/9237 lines (multiplet 1), S I 8693/8694 lines (multiplet 6), S I 
6757 line (multiplet 8), and S I 6046/6052 lines (multiplet 10), separately, which we hereinafter 
referred to as [S/Fe] 92, [S/Fe] 86, [S/Fe] 67, and [S/Fe]6o, respectively. In the case that equivalent- 
width data are available for more than one S I line belonging to the same multiplet, we calculated 
[S/Fe] for each line and adopted their simple mean. 6 

The finally resulting [S/Fe] vs. [Fe/H] relation and the metallicity dependence of the non- 
LTE correction are depicted in figures 5a and b, respectively. Also, the details of these analyses 
(the data of the used equivalent widths and the adopted parameter values, the resulting non- 
LTE abundances or [S/Fe] values with the non-LTE corrections, given for each line/multiplet 
and for each star) are given in electronic table E3 (cf. the footnote in the first page). 

In figure 5a, we can see an interesting trend concerning the behavior of [S/Fe] and the 
importance of the non-LTE effect in metal-poor stars: 

- First, regarding disk stars (—1 < [Fe/H] < 0), the [S/Fe] values of different multiplets almost 
agree with each other (though with a rather large dispersion), showing a gradually increasing 
trend with a decrease of metallicity from [S/Fe] ~ (at [Fe/H] ~ 0) to [S/Fe] ~ 0.3 (at [Fe/H] 
~ -1). 

- However, as the metallicity is further lowered, the discrepancy between [S/Fe] 9 2 and [S/Fe] 8 6 
becomes progressively large, in the sense that [S/Fe] 8 6 continues to rise up to ~ +0.8 (at [Fe/H] 
~ —2.5) while [S/Fe] 92 shows a sign of leveling-off at 0.2-0.3 (or even a sign of downward 
bending at [Fe/H] < -2). 

4.2. 9212/9228/9237 vs. 8693/8694 Discrepancy 

Given the existence of such an apparent discrepancy, we have to contend with a new 
question; i.e., which trend better represents the truth, [S/Fe]g 2 or [S/Fe] 86 ? 

Let us recall here that this difficulty is neither new nor restricted to the present anal- 
ysis including the non-LTE effect. While our non-LTE correction on [S/Fe]g2 (by 0.1-0.3 dex 
acting in the direction of suppressing the S abundance; cf. figure 5b) is surely responsible for 
exaggerating this discrepancy, the non-LTE effect is not the only cause for this discordance. 
As a matter of fact, a quantitative inconsistency between [S/Fe] 92 and [S/Fe] 8 6 has already 
been observed even in pure LTE analyses (see, e.g., figure 4 of Ryde, Lambert 2004; figure 3 
of Takada-Hidai et al. 2005). As mentioned in section 1, however, the recent S abundance 
studies (e.g., Nissen et al. 2004; Ryde, Lambert 2004) tend to preferably invoke [S/Fe] 92 (on 
the assumption of LTE) based on the 9212/9228/9237 lines because of their easy detectability 
even at very low metallicity, while an implicit tendency of neglecting [S/Fe] 8 6 appears to exist, 
presumably because it is based on difficult and less reliable measurements of weak 8693/8694 

6 Regarding multiplet 1 (9212/9228/9237), we did not use the S 1 9228 line in deriving [S/Fe] 92 , because of 
its less reliability due to the blending with H 1 P 9 . 
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lines especially in very metal-poor stars. However, is this really a reasonable attitude? Since 
the situation is rather complicated, let us sort out our thoughts, while reconsidering the validity 
of the so-far adopted assumption. 

First of all, it is worth pointing out that any of the recent arguments suggesting the 
validity of LTE for the S I 9212/9228/9237 lines do not appear to be convincing and should be 
viewed with caution: 

- Nissen et al. (2004) concluded from their analysis on those stars where both multiplet lines 
are measurable that the non-LTE effect should be insignificant because their LTE abundances 
derived from 8693/8694 and 9212/9237 lines turned out to be in good agreement with each 
other (mean difference is 0.03 dex and the standard deviation is 0.08 dex). Actually, we al- 
most confirmed this consequence by our reanalysis of their data. 7 However, the stars they used 
for this check were all in the metallicity range of -1.8 < [Fe/H] < -0.7, and did not include 
any "very" metal-deficient stars at [Fe/H] < —2. This must be a drawback of their argument; 
i.e., the result obtained at —2 < [Fe/H] should not be simply stretched to the region of —3 < 
[Fe/H] < —2 without an actual confirmation. As a matter of fact, as can be seen from fig- 
ure 5a, the A 8 q — A 92 discrepancy begins to be appreciable at the very such region of [Fe/H] 
< —2. As an example, we can refer to the case of HD 140283 ([Fe/H] = —2.42), for which 
they did not perform the A$q vs. A$ 2 check because its EW(8694) was not measured. It can 
be confirmed, however, that A 86 is appreciably larger than A 92 for this very metal-poor star 
even in the LTE assumption, as previously reported in table 3 of Takada-Hidai et al. (2005) 
as (A$ 6 — A 92 )lte = +0.35 (see also the Appendix for a more detailed discussion on the Age 
determination for HD 140283). Accordingly, we consider that it is difficult to justify Nissen et 
al.'s (2004) assertion for the practical validity of LTE . 

- Ryde and Lambert's (2004) statement that "NLTE effects should be small" because "neu- 
tral sulfur atoms represent that main ionization state throughout most of the atmosphere" 
is evidently inadequate, because such a logic can apply only for the case of departures from 
LTE in ionization equilibrium (e.g., the frequently argued case that Fe I lines are affected by 
non-LTE overionization, while Fe II lines are not). As described in subsection 2.1, since the 
non-LTE effect on S I lines is due to the departures from LTE excitation, whether or not S I is 
the dominant ionization stage is irrelevant. 

We would point out, on the contrary, that there is a good reason to suspect the sig- 
nificance of the non-LTE effect, at least for the S I 9212/9228/9237 lines of multiplet 1. The 
transition of this multiplet 1 is 4s 5 S° - 4p 5 P, where the lower 5 S° term (6.52 eV) is "metastable" 
(i.e., its radiative connection to the ground triplet term 3p 4 3 P is only via a very weak forbidden 
transition) and the upper 4p 5 P term is just the lower term of the S I 8693/8694 transition of 
multiplet 6 (4p 5 P - 4d 5 D°) . It should be stressed here that this situation is quite similar to the 

7 For those 19 stars (out of their 34 stars) where both multiplet lines available, we obtained ((As6 — ^92)lte) = 
-0.05(±0.09) and ((A 86 - A, 2 )nlte) = +0.05(±0.08). 
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case of neutral oxygen. Namely, the I 7771/7774/7775 lines of multiplet 1 (3s 5 S° - 3p 5 P) 
and I 6156/6158 lines of multiplet 10 (3p 5 P - 4d 5 D°) just correspond to S I 9212/9228/9237 
and S I 8693/8694, respectively. Then, the formation mechanism of the O I 7771-5 lines may be 
informative for understanding the non-LTE effect of S I 9212/9228/9237. Namely, the pseudo 
two-level-atom nature of the O I 7771-5 line formation (originating from the metastable lower 
level) may approximately also apply to the present case of S I 9212/9228/9237, which may 
suffer an appreciable non-LTE effect depending on the line-strength, such as the case of the 
O I triplet lines at 7771-5 A (see, e.g., Takeda 2003). 8 According to this consideration, we had 
better realize as a starting point that LTE may not be a good assumption for S-abundance de- 
terminations of very metal-poor stars (especially for the triplet lines of multiplet 1 at 9210-9240 
A), while honestly accepting the discrepancy between the LTE abundances of 9212/9228/9237 
and 8693/8694 lines. 

Yet, this does not solve the currently confronted problem, because the non-LTE cor- 
rections derived from our calculations act even in the direction of increasing the discrep- 
ancy. Hence, it is certain that our computed non-LTE corrections (for either or both of 
9212/9228/9237 and 8693/8694 lines) are not adequate, 9 at least for the purpose of appli- 
cations to very metal-poor stars ([Fe/H] < —2), even if they are successful in removing the 
discordance in the LTE abundances of near-solar metallicity F-type stars (cf. section 3). Given 
such a confusing situation, any discussion about the trend (flat? rising?) of the [S/Fe] vs. 
[Fe/H] relation in the very metal-poor regime and its comparison with theoretical chemical 
evolution calculations would not be very meaningful (though several representative theoretical 
trends are overplotted in figure 5a) until the cause of the problem has been clarified. 

4-3. What Is Wrong and What Should Be Done? 

While we can not decide which of the [S/Fe]g2 and [S/Fe]g6 is more reliable, some com- 
ments (if speculative) may be worth presenting here toward a future settlement of this issue: 

- It should be kept in mind that our non-LTE calculations turned out to be successful in rea- 
sonably accomplishing an abundance consistency for the case of F stars of normal metallicity (a 
Per, Polaris, Procyon). Hence, we feel it rather unlikely that some fatally improper modeling 
of the non-LTE calculation (e.g., a serious mis-choice of h, lack of important transitions/levels, 
etc.) is responsible here. Instead, we consider that the cause of the flaw in question should be 

8 One notable difference between the cases of oxygen and sulfur is that, while the LTE ionization equilibrium 
of O i/O II nearly holds due to an efficient charge-exchange reaction leading to a practically LTE population 
for the ground level, such a condition is not expected for the case of S. 

9 More precisely, our non-LTE corrections presented here are not sufficient to reasonably bring the abun- 
dances from different multiplcts into consistency at [Fe/H] < — 2. While this may imply that our non-LTE 
calculations are simply insufficient/incorrect and inapplicable here, an alternative possibility is that some 
non-classical effect (e.g., due to 3D inhomogeneous atmosphere, missing opacity, etc.) other than the non- 
LTE effect may become newly effective at this metallicity range. 
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such that appearing or becoming evident only at the very low metallicity regime. 

- Although we are not qualified to remark on the 3D effect of atmospheric inhomogeneity, 
it seems difficult (at least in a quantitative sense) to invoke this effect in order to remove the 
discrepancy, as long as we see the simulation results presented by Nissen et al. (2004). 

- There might be a possibility of "missing opacity" in our abundance calculation program. 
Namely, if there is some unknown continuum opacity (i.e., not included in the WIDTH9 pro- 
gram we adopted) such that being less sensitive to the metallicity than the H~ opacity, it might 
become significant only in the very metal-poor regime (even if it is overwhelmed by H~ in the 
metal-rich case). Then, the theoretical strengths computed for a given abundance would be 
overestimated, leading to an underestimation of the derived abundances. Such a problem (if 
any exists) might be more probable in the 9210-9240 A region, where a close examination (or 
a photometric matching) of the stellar continuum shape is difficult owing to crowded telluric 
water vapor lines, rather than the much better behaved 8690 A region. In any case, such a con- 
cern should be checked for representative very metal-poor stars by comparing the abundances 
of other elements from this region with those from other regions. 

- We still cannot rule out the possibility that the problem exists in our non-LTE calculations, 
in the sense that the involved errors/flaws may become conspicuous only at the considerably 
low-metal condition. Here, we have a suggestion that might be an important touchstone of our 
calculation. Namely, we recommend to observe the S I 10455/10456/10459 lines of multiplet 3 
for investigating the [S/Fe] behavior at [Fe/H] < —2, which are as strong as 9212/9228/9237 
lines (i.e., usable even for very metal-poor stars) and favorably located in a wavelength region 
almost free from telluric lines. Since these 10455-9 lines belong to the triplet system (4s 3 S° - 
4p 3 P) unlike the quintet-system lines of 9212/9228/9237 and 8693/8694, it may be a validity 
check of our non-LTE modeling to compare the resulting abundances from different systems 
with each other. 

- Though we would refrain from entering the "flat vs. rising" controversy on the [S/Fe] trend, 
it should be kept in mind that our understanding of the formation S I 9212/9228/9237 lines is 
still incomplete, which thus should not be trusted too much. Meanwhile, we should not forget 
the merit of the weaker and deep- forming 8693/8694 lines (as far as measurements can be ac- 
curately done), to which the classical line-formation theory may be applied more safely. In this 
sense, the "rising" tendency suggested from [S/Fe]s6 (Israelian, Rebolo 2001; Takada-Hidai et 
al. 2002) might deserve more attention. 

- In our opinion, what we need is to carry out much more "reliable" observations of weak S I 
8693/8694 lines for as many very metal-poor stars ([Fe/H] < —2) as possible, since such data 
are evidently insufficient. Admittedly, the difficulty of such observations would progressively 
increase as we go into such a low- metal region, and spectra of very-high S/N would be required. 
It should be here mentioned that very metal- deficient halo stars of comparatively higher T e g 
are suitable for this purpose, since the strengths of such high-excitation S I lines increase with 
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T e ff. Hence, as a possibility, it may be promising to pay attention to F-type horizontal-branch 
stars in very metal-poor globular clusters (e.g., M92; [Fe/H] = —2.3). 

5. Behavior of Zinc in Disk/Halo Stars 

Similarly to the case of sulfur, we also carried out an extensive non-LTE reanalysis 
of the published equivalent-width of Zn I 4722/4810 and 6362 lines taken from the following 
papers: 10 Sneden and Crocker (1988), Sneden, Gratton, and Crocker (1991), Beveridge and 
Sneden (1994), Prochaska et al. (2000), Nissen et al. (2004), Cayrel et al. (2004), Honda et al. 
(2004), and Chen, Nissen, and Zhao (2004). In a similar way as described in subsection 4.1, 
non- LTE Zn abundances (A^ TE ) were determined from these EW data of 4722/4810/6362 11 
lines along with the model atmospheres corresponding to the atmospheric parameters taken 
from the same papers of the EW source. The [Zn/Fe] ratio was derived as 

[Zn/Fe] = « LTE - 4.60) -[Fe/H], (2) 

where 4.60 is the solar zinc abundance (Anders, Grevesse 1989). Treating the Zn I 4722/4810 
lines (multiplet 2) and the 6362 line (multiplet 6) separately, we derived [Zn/Fe] 4722/4810 and 
[Zn/Fe] 6362- In the case where both of the 4722 and 4810 lines are available, we adopted a 
simple average of the two to obtain [Zn/Fe] 472 2/48io- The finally resulting [Zn/Fe] vs. [Fe/H] 
relation and the metallicity-dependence of the non-LTE correction are depicted in figures 6a 
and b, respectively. As in the case of S, the details of these Zn reanalyses are given in electronic 
table E4 (cf. the footnote in the first page). 

By inspecting figures 6a and b, we can see the following characteristics: 

- Generally speaking, the non-LTE corrections in zinc abundance determinations from Zn I 
4722/4810 and 6362 lines are comparatively insignificant 0.1 dex in most cases). Especially, 
LTE is a practically valid approximation for the Zn I 6362 line, which is appropriate for inves- 
tigating the Zn abundances of disk stars as recently used by Chen, Nissen, and Zhao (2004). 
The corrections are mostly positive (i.e., the non-LTE effect tends to act in the direction of 
weakening the line) except for some cases of the 4722/4810 lines, and have a slightly increasing 
trend with a lowering of the metallicity, which gives a mild influence on the [Zn/Fe] vs. [Fe/H] 
relation as described below. 

- The traditional belief that "Zn almost scales with Fe over a wide range of metallicity (i.e., 
[Zn/Fe] ~ 0)" does not adequately represent the truth from a strict point of view. Namely, the 

10 We remark that there are other important recent papers apart from these eight concerning Zn abundance 
analyses for a number of stars, such as Primas et al. (2000), Mishcnina et al. (2002), Reddy et al. 
(2003), Bcnsby, Feltzing, and Lundstrom (2003), Allcndc Prieto et al. (2004), and Ecuvillon et al. (2004). 
Unfortunately, since these authors presented only the abundance results and did not publish observed EW 
values for each star, we could not include their data in our non-LTE reanalyses. 

11 In deriving the abundance from the Zn 1 6362 line, we did not take into account the broad wing of the Ca I 
autoionization line, since it is practically negligible (cf. Chen et al. 2004). 
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[Zn/Fe] values for disk stars gradually increase from [Zn/Fe] ~ (at [Fe/H] ~ 0) to [Zn/Fe] ~ 
0.2 (at [Fe/H] ~ —1). While it appears that a kind of weak discontinuity exists at [Fe/H] ~ — 1, 
[Zn/Fe] exhibits a pseudo-plateau at ~ 0.2 (or a slightly increasing trend with a very gentle 
slope) over the region of —2 < [Fe/H] < —1. Then, below [Fe/H] ~ —2, [Zn/Fe] begins to show 
a manifest rise toward an extremely low metallicity regime, which continues down to [Fe/H] 
~ —4 where [Zn/Fe] attains even a value of ~ +1. In short, we have simply confirmed the 
recently reported results, such as the slight increase of [Zn/Fe] for disk stars with a decrease of 
[Fe/H] (Reddy et al. 2003; Chen et al. 2004) and the clear evidence of markedly rising [Zn/Fe] 
at [Fe/H] < —2.5 (Primas et al. 2000; Cayrel et al. 2004), while these trends have become even 
more pronounced in our reanalyses because of the (positive) non-LTE corrections gradually 
increasing with a decrease of metallicity. 

- It appears that none of the theoretical predictions from representative Galactic chemical 
evolution calculations (such those overplotted in figure 6a) adequately explain such a trend of 
[Zn/Fe] with a run of [Fe/H], though (somewhat interestingly) only the calculation by Goswami 
and Prantzos (2000) under the assumption of "constant yield fixed at the solar metallicity" (i.e., 
labeled "GP(const)"; an unrealistic case only for an illustration purpose to show the signifi- 
cance of the metallicity-dependence of the yield) appears to roughly show a kind of similar 
tendency (though only qualitatively, not quantitatively). New theoretical calculations that 
can adequately reproduce such an observationally established [Zn/Fe] vs. [Fe/H] relation are 
awaited. 

6. Concluding Summary 

We carried out statistical-equilibrium calculations on neutral sulfur and zinc for an 
extensive grid of models over a wide range of parameters (T e g- from 4500 K to 7000 K, logg 
from 1.0 to 5.0, and [Fe/H] from —4.0 to +0.5), in order to study how the non-LTE effect is 
important in deriving the photospheric S and Zn abundances of F, G, and K stars from various 
S I and Zn I lines, where a particular emphasis was placed on S I 8693/8694 lines (multiplet 6), 
S I 9212/9228/9237 lines (multiplet 1), Zn I 4722/4810 lines (multiplet 2), and Zn I 6362 line 
(multiplet 6). We then constructed an extensive grid of non-LTE abundance corrections, which 
should be applied to the conventionally derived LTE abundance to obtain the corresponding 
non-LTE abundance. 

Roughly speaking, the non-LTE effect on S I lines acts in the direction of strengthening 
the lines (i.e., negative non-LTE correction) owing to the overpopulation of the lower level 
coupled with the dilution of the line source function. On the other hand, the line source 
function for the Zn I lines tends to be superthermal (Sl > B), while the lower level shows either 
non-LTE overpopulation or underpopulation depending on the cases; as a result, Zn I lines 
tend to be weakened (i.e., positive correction) by the non-LTE effect in many cases, though the 
opposite can occur in some cases. Quantitatively, the extents of the non-LTE corrections for 
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the S I 9212/9228/9237 lines are considerably larger and more important than those for the S I 
8693/8694 lines, while the non-LTE effect for the Zn I lines is generally of minor importance. 

Taking account of the fact that the non-LTE effect tends to become larger with an 
increase/lowering of T e s /log g, we performed abundance analyses of sulfur and zinc for rep- 
resentative F supergiants/subgiant (a Per, Polaris, Procyon) along with the Sun, in order to 
check the validity of our non-LTE calculations by examining whether a consistency can be 
achieved between the abundances derived from different lines. For this purpose, we used the 
high-dispersion echelle spectra obtained with the GAOES spectrograph at Gunma Astronomical 
Observatory. It was confirmed that the large discrepancies seen in the LTE S abundances of 
these F stars could be successfully removed by our non-LTE corrections, while the non-LTE 
corrections for the Zn I lines were too small to be useful for such a check. 

Finally, extensive non-LTE reanalyses of published equivalent-width data of the S I and 
Zn I lines were carried out, in order to investigate the behavior of [S/Fe] or [Zn/Fe] with a 
change of [Fe/H] in Galactic disk/halo stars. The following conclusions were reached from this 
restudy: 

- We encountered a serious difficulty in the [S/Fe] vs. [Fe/H] relation at the very metal-poor 
region of [Fe/H] < —2, in the sense that the [S/Fe] values derived from S I 8693/8694 tend to 
rise progressively with a decrease of metallicity while those from S I 9212/9228/9237 show a 
flat plateau (or a sign of slight downward bending); the discordance amounts up to ~ 0.5 dex. 
Though such a trend has actually been reported even in the framework of LTE, this discrep- 
ancy has been exaggerated due to the application of non-LTE corrections. At present, little 
can be said about which represents the truth. Yet, we have a rather conservative feeling that 
the solution based on the weaker and deep-forming 8693/8694 lines (ever-rising [S/Fe]) might 
deserve more attention (in contrast to the arguments of recent studies), since our understanding 
of 9212/9228/9237 lines is still incomplete, with which unclarified problems might be involved 
(e.g., flaws in non-LTE calculations, 3D effect, missing opacity, etc.). 

- Inspecting the [Zn/Fe] vs. [Fe/H] relation resulting from our reanalysis, we almost con- 
firmed the recently reported tendency, such as the gradual increase of [Zn/Fe] from [Zn/Fe] 

~ (at [Fe/H] ~ 0) to [Zn/Fe] ~ 0.2 (at [Fe/H] 1), nearly constant [Zn/Fe] (i.e., without 

particular trend) over the region of —2 < [Fe/H] < — 1> an d a beginning of rise at [Fe/H] < —2 
continuing toward an extremely low-metallicity regime. Actually, these trends appear even 
more exaggerated in our reanalysis results, because the (positive) non-LTE corrections acting 
on the Zn abundances are metallicity-dependent (i.e., progressively increasing with a decrease 
of metallicity). 

Appendix. Sulfur 8694.63 Line in the Spectrum of HD 140283 

Since the S abundance determination of very metal-poor stars ([Fe/H] < —2) based on 
the S I 8694 line is very difficult because of its considerable weakness, the rare case of HD 140283 
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is especially important, for which one may manage to detect and measure this line as recently 
carried out by Takada-Hidai et al. (2005). Since only the figure in a magnified scale that they 
presented for demonstrating the rather delicate S I 8694 detection (cf. their figure 2d) is not 
necessarily sufficient for the reader to judge its reliability, we show here some supplementary 
and more informative figures. 

As Takada-Hidai et al. (2005) did, we invoked the ESO/UVES spectrum of HD 140283 
in the published high-dispersion stellar spectral library, "A Library of High-Resolution Spectra 
of Stars across the Hertzsprung- Russell Diagram" (Bagnulo et al. 2003). Regarding the cal- 
culation of theoretical spectra to be compared with observations, we adopted the atmospheric 
parameters of T c g = 5960 K, logg = 3.69, v t = 1.5 km s _1 , and [Fe/H] = —2.42, which were 
taken from Nissen et al. (2004). Since the non-LTE effect is practically negligible for such a 
weak S I 8694 line (cf. electronic table E3), we assumed LTE in the spectrum synthesis. A 
comparison of the observed spectrum with three computed spectra corresponding to [S/Fe] = 
0.0, +0.5, and +1.0 is displayed in figures 7a (wide view) and b (magnified view). The theo- 
retical spectra are convolved with a Gaussian broadening function, which was so chosen as to 
accomplish the best fit for the conspicuously seen Fe I 8688.62 line (cf. figure 7a). 

Inspecting these figures, we can see that a weak (but recognizable) dip with a depth of 
~ 1% surely exists at the position of S I 8694.63. Admittedly, we cannot rule out a possibility 
that this is nothing but a fluctuation of fringe patterns. However, since the S/N ratio of this 
spectrum is estimated to be ~ 500 (a ~ 0.2%) from the line-free 8690-8692 A region, the 
possibility of such a large fluctuation (amounting to ~ 5cr) is not considered to be very likely. 
Then, on the standpoint that this identification is real, we may state that [S/Fe] should be near 
to ~ +1, since this line would not be visible if [S/Fe] < +0.5 (figure 7b). 
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Table 1. Atomic data of the relevant S I and Zn I lines. 



Species 


Line 
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8.86 
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( — 7.57) 


S I 
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3 
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—0.430 
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3 
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6.860 


+0.040 
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— 5.21 


( — 7.57) 


S I 
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8693.931 
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—0.510 


7.62 


—4.41 


( — 7.30) 


S r 
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6 


4p P 3 - 4d D 4 
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+0.080 


7.62 


—4.41 


( — 7.30) 


S i 


6757 


8 


A 5 "T> £ ] 5 T~\ O 

4p F3 - 5d D 9 
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All data are were taken from Kurucz and Bell's (1995) compilation as far as available. RMT is the multiplet number given by the Revised 
Multiplet Table (Moore 1959). Gammar is the radiation damping constant, log7 ra d- Gammas is the Stark damping width per electron 
density at 10 4 K, log(7 c /iV c ). Gammaw is the van der Waals damping width per hydrogen density at 10 4 K, log(7 w /iVH). Note that the 
values in parentheses are the default damping parameters computed within the Kurucz's WIDTH program (cf. Lcusin, Topil'skaya 1987), 
because of being unavailable in Kurucz and Bell (1995). The meanings of other columns arc self-explanatory. 
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Table 2. Dependence of the non-LTE effect of S I 8694 and 9212 lines on the H I collision and photoionization cross section. 
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Columns 1—6 are self-explanatory (the units of T c f f , g, and £ arc K, cm s , and km s , respectively). While W in the 7th column is the LTE equivalent 

width (calculated for the atmospheric parameters and the input abundance given in columns 1—5) , the Ws in columns 8—12 and As in columns 13-17 arc the 
non-LTE equivalent width (in m A) and the non-LTE abundance corrections (in dex) , respectively, where the suffixes ( + 1) 0, — 1, — 2, — 3, and — 4) denote the 
corresponding values of h (the logarithm of the H I collision correction factor applied to the classical formula). The values for the finally adopted h = case 
arc highlighted by boldface characters. In the 18th and 19th columns arc given the variation of Aq caused by chaging the photoionization cross section (for all 
levels) by a factor of 1/10 [<5(A_) = A (a atd ) - A (a atd X 0.1)] and 10 [S(A + ) EE An^a^) - A (e« 3td X 10)]. 
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Table 3. Dependence of the non-LTE effect of Zn I 4810 and 6362 lines on the H I collision and photoionization cross section. 
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Sec the note in tabic 2 for a detailed description of the presented data. The non-LTE Ws and As could not be successfully calculated for the case of (T Q ff = 4500 
K, logg = 2.0, and [Fc/H] = -1.0) because of instability problems. 
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Table 4. Non-LTE analysis of S I and Zn I lines for a Per, Polaris, Procyon, and the Sun. 



Line 


wi 


logT 1 " 


.NLTE 
Yi 


a LTE 
A 


A* 


<5 vd„+ 


a Per 


(T eff = 6250, log 5 = 


- 0.90, [A J 


- 0.0, £ - 


- 4.5) 




o t o c n o i A 

b I 8693+4 


264.7 


-0.77 


6.96 


7.21 


— 0.25 


0.00 


S I 9212 


521.4 


-2.36 


7 1/1 

i . r^ 


8.24 


— 1 .10 


— 0.01 


C t nno7 

b I 9z3 ( 


455.7 


-2.18 


7.22 


8.31 


— 1.09 


0.00 




80.1 


-0.62 


4.10 


4.01 


+0.09 


0.00 


Zn I 4810 


87.8 


-0.68 


4.01 


3.91 


+0.10 


0.00 


Zn I 6362 


35. 8 11 


-0.42 11 


4.48 11 


4.60 11 


-0.12 11 


0.00 


Polaris 


(T eff = 6000, logg; 


— 1.50, [X] ■■ 


— 0.0, £ ■ 


— 5.0) 




S I 8693 


68.0 


-0.50 


7.04 


7.19 


—0.15 


0.00 


S I 8694 


128.2 


-0.70 


6.87 


7. 10 


— 0.23 


0.00 




458.6 


-2.04 


D.o I 


l .oil 


— U.J.j 


n nn 
u.uu 


b I 9237 


352.9 


-1.70 


6.82 


7.50 


— 0.68 


0.00 




99.8 


-0.76 


4. 13 


4.U< 


+0.06 


0.00 


Zn I 4810 


110.5 


-0.83 


4.05 


3.99 


+0.06 


0.00 


Zn I 6362 


28.8 


-0.37 


4.34 


4.41 


-0.07 


0.00 


Procyon 


(Teff = 


6600, logg 


a nn r vl 

— 4.0U, [X\ 


— 0.0, £ 


— 2.0) 




S I 8693 


28.8 


-0.51 


7.18 


7.22 


— 0.04 


— 0.01 


S I 8694 


72.2 


-0.78 


7.24 


7.31 


— 0.07 


— 0.04 


S I 9212 


207.7 


-1.93 


7 08 


7 47 


q 


Q 


S I 9237 


161.8 


-1.74 


7. 16 


7.54 


— 0.38 


— 0.07 


Zn I 4722 


65.2 


-0.94 


4.46 


4.49 


— 0.03 


— 0.01 


Zn I 4810 


71.0 


-1.03 


4.38 


4.44 


— 0.06 


— 0.01 


Zn I 6362 


18.8 


-0.40 


4.46 


4.52 


— 0.06 


0.00 


Sun (T eff = 5780, logg = 


4.44, [X] — 


0.0, ^ — 


1.0) 




S I 8693 


10.6 


-0.36 


7.18 


7.19 


-0.01 


-0.03 


S I 8694 


28.5 


-0.45 


7.17 


7.18 


-0.01 


-0.06 


S I 9228 


95.1 


-1.07 


7.06 


7.21 


— 0.15 


— 0.15 


S I 9237 


97.1 


-0.95 


7. 14 


7.25 


— 0.11 


— 0. 16 


S I 10455 


111.9 


-0.82 


f . 11 


7.20 


— 0.09 


— 0. 16 


S I 10456 


55.3 


-0.58 


7.14 


7.19 


-0.05 


-0.08 


S I 10459 


88.2 


-0.73 


7.09 


7.16 


-0.07 


-0.10 


Zn I 4722 


67.4 


-1.23 


4.56 


4.61 


-0.05 


-0.13 


Zn I 4810 


71.6 


-1.28 


4.49 


4.54 


-0.05 


-0.16 


Zn I 6362 


20.5 


-0.50 


4.53 


4.53 


0.00 


-0.03 



* Equivalent width in units of mA. 

^ Mean line-formation depth in terms of the standard continuum optical depth at 5000 A, which was calculated (for the non-LTE case) in 
the same manner as described in Takeda and Takada-Hidai (1994). 
x Non-LTE correction defined as A ~ J 4 NLTE - A LTE . 

§ Variation of ^4 NLTE caused by increasing the van dcr Waals damping width (for which we adopted the default treatment of the WIDTH9 
program, equivalent to the classical Unsold's approximation; cf. table 1) by a factor of 2.5. 

" These values are less reliable and should be viewed with caution because of the difficulty in measuring the equivalent width (cf. figure 
3c). 
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Fig. 1. Ratio of the S I line source function to the local Planck function (B) and the NLTE-to-LTE 
line-center opacity ratio as functions of the standard continuum optical depth at 5000 A computed for 
models of T c g — 4500 K, 5500 K, and 6500 K. The green lines and blue lines correspond to S'l/B and 
Zq' lte //q TE , respectively. The solid lines show the results for the 4s 5 S° - 4p 5 P transition of multiplct 1 
(corresponding to S I 9212/9228/9237 lines), while those for the 4p 5 P - 4d 5 D° transition of multiplet 
6 (corresponding to S I 8693/8694 lines) are depicted" by dashed lines. In each case, the results for two 
different gravity atmospheres are given: The thick lines are for logg = 4 and the thin lines are for logg = 2. 
Note also that the curves are vertically offset by an amount of 0.5 dex relative to those of the adjacent 
metallicity ones. 
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Fig. 2. Ratio of the Zn 1 line source function (Sl) to the local Planck function (B) and the NLTE-to-LTE 
line-center opacity ratio as functions of the standard continuum optical depth at 5000 A computed for 
models of T c g — 4500 K, 5500 K, and 6500 K. The green lines and blue lines correspond to S^/B and 
Iq LTE / Iq TE , respectively. The solid lines show the results for the 4p 3 P° - 5s 3 S transition of multiplct 
2 (corresponding to Zn 1 4722/4810 lines), while those for the transition 4p 1 P° - 4d 1 D of multiplct 6 
(corresponding to Zn 1 6362 line) are depicted by d^ihed lines. In each case, the results for two different 
gravity atmospheres are given: The thick lines are for \ogg = 4 and the thin lines are for logg = 2. Note also 
that the curves are vertically offset by an amount of 0.5 dex relative to those of the adjacent metallicity 
ones. 




Wavelength (A) 

Fig. 3. Spectra of three bright F stars observed by using the GAOES spectrograph at Gunma 
Astrophysical Observatory at five wavelength regions (a — Zn I 4722 line region, b — Zn I 4810 line re- 
gion, c— Zn I 6362 line region, d— S I 8693/8694 lines region, e— S I 9212/9228/9237 lines region), on 
which the equivalent widths of these S I and Zn I lines (their positions are indicated by downward ar- 
rows) were measured for an adequacy check of our non-LTE calculations. The KPNO solar flux spectra 
of Kurucz et al. (1984) are also shown for comparison. The spectra are placed according to the order of 
a Per, Polaris, Procyon, and the Sun from top to bottom, each being vertically offset by an appropriate 
constant (0.2, 0.3, 0.075, 0.2, and 0.4 for panels a, b, c, d, and e, respectively) relative to those of the 
adjacent metallicity ones. Note that, in panel (e) of 9210-9240 A region, numerous telluric lines due to 
H2O have been removed for the GAOES spectra by dividing them by the spectrum of a rapid rotator (7 
Cas), unlike the KPNO solar spectrum where those telluric lines arc conspicuously observed. 
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9236 9237 9238 9239 

o 

Wavelength (A) 

Fig. 4. Observed profiles (open circles; the data are the same as figure 3) of the S I line at 9237.538 A for a 
Per, Polaris, Procyon, and the Sun (from top to bottom; each spectrum is vertically shifted by 0.5 relative 
to the adjacent one), fitted with the theoretically calculated profiles (solid lines). The theoretical profiles 
were computed with the non-LTE abundances (^4 NLTE ) derived from the equivalent- width analysis along 
with the atmospheric parameters presented in table 4, and then convolved with Gaussian broadening 
functions appropriately chosen so as to make the best fit. (The continuum levels and the wavelength 
scales of the observed spectra have also been so adequately adjusted as to accomplish the best match.) In 
addition to the non-LTE profiles depicted in thick lines, the corresponding LTE profiles (computed also 
with ^nltE) are also shown in thin lines. For the daQe of the Sun, since the overlapping wings of telluric 
lines are not included in our spectral synthesis, the fit does not appear to be satisfactorily good. The weak 
absorption feature at A ~ 9238A recognized in the theoretical solar spectrum is due to the Si I 9238.04 
line, which is also blended with a strong telluric line. 




Fig. 5. (a) [S/Fe] vs. [Fe/H] relation resulting from our non-LTE reanalysis of the published equiva- 
lent-width data of S I lines taken from various literature. Open circles — results from S I 9212/9237 lines 
of multiplet 1; filled circles — results from S I 8693/8694 lines of multiplet 6; Greek crosses (+) — results 
from S I 6756 line of multiplet 8; St. Andrew's crosses (x) — results from S I 6052 line of multiplet 10. Note 
that the larger symbol corresponds to low-gravity giant stars (logg < 3) and the smaller symbol corresponds 
to high-gravity dwarf /subgiant stars (logg > 3). The representative theoretical predictions are depicted 
by lines: Dash-dotted line (S98)— taken from figure 12 of Samland (1998); solid lines (TWW 1/2 x, lx, 
2x) — taken from figure 22 of Timmes, Woosley, and Weaver (1995) corresponding to three choices of 
the adjustment factor (1/2, 1, and 2) for the Fc yield from massive stars by which the standard Woosley 
and Weaver's (1995) yield is to be multiplied; dashed/dotted lines — taken from figure 7 of Goswami and 
Prantzos (2000) for the two cases of S yield, i.e., the dotted line is for the metallicity-independent yield 
[GP (const)] and the dashed line is for the realistic mctallicity-dependent yield [GP(var)]. (b) The corre- 
sponding non-LTE corrections used for deriving the [S/Fe] values shown in panel (a), plotted as functions 
of [Fe/H]. 
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Fig. 6. (a) [Zn/Fe] vs. [Fc/H] relation resulting from our non-LTE reanalysis of the published equiva- 
lent-width data of Zn I lines taken from various literature. Open circles — results from Zn I 4722/4780 
lines of multiplet 2; filled circles — results from Zn I 6362 line of multiplet 6; Note that the larger sym- 
bol corresponds to low-gravity giant stars (logg < 3) and the smaller symbol corresponds to high-gravity 
dwarf /subgiant stars (log<? > 3). The representative theoretical predictions are depicted by lines: Solid 
lines (TWW 1/2 x, lx, 2x) — taken from figure 35 of Timmes, Woosley, and Weaver (1995) corresponding 
to three choices of the adjustment factor (1/2, 1, and 2) for the Fe yield from massive stars by which 
the standard Woosley and Weaver's (1995) yield is to be multiplied; dashed/dotted lines — taken from 
figure 7 of Goswami and Prantzos (2000) for the two cases of Zn yield, i.e., the dotted line is for the 
metallicity- independent yield [GP (const)] and the dashed line is for the realistic metallicity-dependent 
yield [GP(var)]. (b) The corresponding non-LTE corrections used for deriving the [Zn/Fe] values shown 
in panel (a), plotted as functions of [Fe/H]. 
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Fig. 7. Open circles: Spectrum of HD 140283 observed with ESO/UVES, which was taken from the 
spectral database of Paranal Observatory Project (Bagnulo et al. 2003). Solid lines: Theoretical spectra 
computed with the atmospheric parameters (T e g, logg, v t , [Fc/H]) of (5690 K, 3.69, 1.5 km s -1 , —2.42) for 
three sulfur abundances of [S/Fc] = 0.0, +0.5, and +1.0, where the calculation was done in LTE, because 
the non-LTE effect is negligibly small for such very weak S I 8693-4 lines at this parameter range (cf. 
electronic table E3). The computed spectra were convolved with a Gaussian function, adequately chosen 



